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After constructing all the tetraquark interpolating currents with J^'~^ = 1""^, l'^"'' and 1^ 
in a systematic way, we investigate the two-point correlation functions to extract the masses of the 
charmonium-like states with QCD sum rule. For the 1^^ qcqc charmonium-like state, mx = 4.6 ~ 
4.7 GeV, which implies a possible tetraquark interpretation for the state Y'(4660). The masses for 
both the l"""^ qcqc and scsc charmonium-like states are around 4.0 ~ 4.2 GeV, which are slightly 
' above the mass of X(3872). For the 1 and 1^ qcqc charmonium-like states, the extracted masses 

are around 4.5 ~ 4.7 GeV and 4.0 ~ 4.2 GeV respectively. As a byproduct, the bottomonium-like 
, states are also studied. We also discuss the possible decay modes and experimental search of the 

^ ■ charmonium-like states. 
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I. INTRODUCTION 



In the past several years, many unexpected charmonium-like states have been discovered at B-factories, some of 
which lie above the open charm threshold and decay into final states that contain a cc pair. Some of them do not 
fit in the conventional quark model easily and are considered as the candidates of the exotic states beyond the quark 
model, such as the molecular states, tetraquark states, the charmonium hybrid mesons, baryonium states and so on. 
For experimental reviews of these new states, one can consult Refs. [iHl]. 

The underlying structure of these new states inspired the extensive study of the hadron spectroscopy. X(3872) is 
, the best studied charmonium-like state since its discovery by the Belle Collaboration 0. Although the analysis of 
■ angular distributions favors the assignment J^'^ = 1++ [3,3, the 2 ^ possibility is not ruled out Q. The mass and 
decay mode of AT (3872) are very different from that of the 2^ Pi cc state. Up to now, the possible interpretations 
of AT (3872) include the molecular state [l0l - [l6| . tetraquark state [l7l - [l9| . cusp [l^l and hybrid charmonium [2l[. In 
Ref. the authors studied the J^*-^ = 1++ state using a tetraquark current in the framework of the QCD sum rule 
approach. 

The initial state radiation (ISR) process played an important role in the search of the 1 charmonium-like states 
at B-factorics. BaBar Collaboration first observed y(4260) in the e+e~ -> 'yisRJ/4'T^'^T^~ process [l^l, which was 
' confirmed by Belle Collaboration [2^. Then, BaBar studied the ^isRip{2S)Tr~^n~ channel and observed a broad 
^ I enhancement around 4.32 GeV. Using the same technique. Belle observed two distinct resonances F(4360) and F(4660) 



in the V'(2S')7r+7r~ mass distribution [2^. The masses of these new charmonium-like states are higher than the open 
charm threshold. However, the Y — >■ decay modes have not been observed yet [2^, which are predicted 

to be the dominant decay modes of the charmonium above the open charm threshold in the potential model. In 
Refs. [l^H^l, the authors studied the 1 charmonium-like Y mesons u sing the QCD sum rule approach. Maiani 
et al. tried to assign F(4260) as the scsc tetraquark in a P-wave state [231 ■ ^(4260) was also interpreted as the 
interesting charmonium hybrid state [2^-[3l|. F(4660) was considered as a V'(2S')/o(980) bound state 

The exotic state with J^*^ = 1 can not be a qq state in the simple quark model. Neither can the 1 ^ state be 
formed by two gluons due to the Landau- Yang selection rule [33l |33 | . States with such a quantum number a re g ood 
candidates of the hybrid meson, which has been studied with the MIT Bag model [HI, [s^ , the flux tube model [13 Issj 
and the QCD sum rule formalism [3^|4^. Recently there have been some efforts on the 1 ^ charmonium-like exotic 
states. For example, the structure of Ar(4350) was studied using a D*D*q current with J^'~'' = 1 4l|. Moreover, 
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the newly observed state y(4140) was argued as a 1 exotic charmonium hybrid state [42| . 

In Ref. [1^, Nielsen et al. reviewed the charmonium-like states from the perspective of the QCD sum rule ap- 
proach. They studied new resonances such as X(3872), Z+(4430) and F(4660) etc using one single molecular-type 
or tetraquark-type interpolating current. In this paper, we first construct all the local charmonium-like tetraquark 
currents with J^'-^ = 1 ^, 1 , 1++ and 1"' in a systematic way. With these independent currents, we study the 
two-point correlation functions and extract the masses of the possible 1 1^^, 1^+, 1"' states. We study both the 
qQqQ and sQsQ systems where Q ~ c,b. 

The paper is organized as follows. In Sec. II, wc construct the tetraquark interpolating currents with J^'~' — 
1 ^, 1 , 1++ and 1^ using the diquark and antidiquark fields. In Sec. Ill, we calculate the correlation functions 
and spectral densities of the interpolating currents and collect them in the Appendix. We perform the numerical 
analysis and extract the masses in Sec. IV and discuss the possible decay modes and experimental search of the 
charmonium-like states in Sec. V. The last section is a brief summary. 



II. TETRAQUARK INTERPOLATING CURRENTS 



In this section, we construct the diquark- antidiquark type of currents using the technique developed in our previous 
works [43 - l46j . Considering the Lorentz structures, there are five independent diquark fields (or anti-diquark fields): 
q^Cj^Cb, q^Ccb, qaCjf^j^ct, q'^C'^f^Cb and q^Ca^^Cb, where a, b are color indices and q denotes an up or down quark. 
We will also take account of q^Ca^u^^Cb although it is equivalent to q^Ca^j^yCb- In fact, they have different parities. 
Using these diquarks and anti-diquarks as the basis, we can compose a six-order matrix O. The elements of O are 
the tetraquark operators without color structures. We show the spins and parities of the matrix elements with J < 1 
in Tabled 



Operators 




qaibCcl' 


qaCcl' 


qait^ibCcl' 


qait^Ccl' 








JP 


0+ 


Q- 


l- 


1+ 




1+ 




0+ 


0+ 


Q- 




1+ 






qlCcb 


0- 


0- 


0+ 


1+ 


1- 






qICjfj,"f5Cb 


1- 




1+ 


0+ 




1+ 




qlC-yf.Cb 


1+ 


1+ 






0+ 




1+ 


qJCafj^^Cb 








1+ 




0+ 




qaCa^^^jsCb 


1+ 






1- 


1+ 







TABLE I: The spins and parities of the elements of the matrix O. 



We do not use Oee and Oes in the construction of the currents since Oee is equivalent to O55 and Ogs equivalent 
to O56 . One notes that under the charge-conjugation transformation: 

COC-i = O^, (1) 

where is the transpose of O. With this relation, we can define the symmetric matrix S and antisymmetric matrix 
A: 

S = + 0^,A = 0-0^ (2) 

The elements of these two matrices are the tetraquark operators that have definite C-parity: they have even and 
odd C-parities for the elements of 5* and A, respectively. An = indicates that the J^*^ = 0"' tetraquark currents 
without derivatives do not exist, which has been proven in Ref. [45| . 

The diquark and antidiquark should have the same color symmetries to compose a color singlet tetraquark current. 
So the color structure of the tetraquark is either 6 Cg) 6 or 3 (g) 3, which is denoted by 6 and 3 respectively. Finally, we 
can obtain the tetraquark interpolating currents with J^*^ = 1 ^ , 1^^, 1++ and 1^ from the matrices S and A: 
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• The interpolating currents with J^'~' = 1 and 1 are: 

-^iM = Sly,{Al^) = qaCj5Cb{qalf,75Ccl + qblt^lbCcl) ± qaCj^j^CbiqalbCcl + qblbCcl) , 

J2t^ = 515(^15) = qlC-i^'cbiqaf^^uCcl - qbUy^yCcl) ± qlCa^^Cbiqal^Ccl - qbi^Ccl) , 

•^3m = SfsiAl^) = qaCj5Cb{qalt.l5Ccl - qblf^lbCcl) ± qaC-ff,-f5Cb{qal5Ccl - qb-f5Ccl) , 

Jip- = -S'lsl^ts) = qICYcb{qaa^„Ccl + qba^^Cc^) ± qlC(Ji,^Cb{qaYCcl + qb-f^Ccl) , 

J5ti = S24{A^4) = qaCcb{qa7fj.Cd[ + qbj^Ccl)±q'^Cj^Cb{qaCcl + qbCcl), (3) 

•^6m SIq{A%) = qlC-i^-izCbiqa^^^LvlbCcl + qb(^ t,i,l^Cc^) ± grCcr^j,75Cf,(ga7''75(:7c^ + qbl^lbCc^) , 

J?,. = -524(^24) = qlCcbiqalt^Ccl - qbJf^Ccl) ± qJCj^,Cb{qaCcl - ^feCc^) , 

•^Sm = 'S'|6(-436) = qlCj^-fsCbiqaCTp^^JbCcl - gh(T^^.75Cc^) ± qlCa^j,y-i^Cb{qaYlbCcl ~ qbY-fsCcl) . 
where "S"' and "+" correspond to J^'-' = 1 ^, "A" and "— " correspond to J^'-' — 1 . 

• The interpolating currents with J^'-' = 1++ and 1"' are: 

•^iM = = qICcb{qalt,l5Ccl + qblt^lbCcl) ± qlC-if,-i5Cb{qaGcl + qbCcl) , 

•^2^ = S23{Als) = qaCcb{qalt.l5Ccl - qblf^lbCcl) ± qlC-if,-i5Cb(qaCcl - qbCcl) , 

Jsfi = 'S'i4(^i4) = qlC-izCb{qaltiCcl + qbjf,Ccl) ± qlC-iy,Cb{qalbCcl + gfc75Cc^) , 

>/4m = 'S'?4(^i4) = qlC-i^Cb{qaltiCcl - qbjf,Ccl) ± qlC-iy,Cb{qalbCcl - qblbCcl) , (4) 

-'"sm = 'S'46(^46) = qlCl" Cb{qaO t^vlbCcl + gh(T^^75Cc^) ± qlCa^,^-i^Cb{qaYCcl + qbj^Ccl) , 

-'^eM = '5'|6(^46) = qlCl^CbiqafTt^vlbCcl - qbC^t^vlhCcl) ± qlCa^,^-i^Cb{qaYCcl - qbj^Ccl) , 

-^Ta. = '535(^35) = qlCl''lbCb{qaOy,vCcl + qbUy^yCcl) ± qlCuy^yCbiqal^lbCcl + qbYlbCcl) , 

-^Sm = 'S'|5(A35) = q^CYj5Cb{qa<y,j.uCcl - qbUy^yCcl) ± qlCa^^yCbiqal^lbCcl - qbYlbCcl) . 

where "5" and "+" correspond to J^'-^ = 1++, "A" and "— " correspond to J^'-' = 1^ . 

In order to have definite isospin and G-parity, all the currents in Eqs. ©-(U) should contain (ucuc + dcdc). However, 
we do not differentiate the up and down quarks in our analysis due to the isospin symmetry and denote them by q. 



III. SPECTRAL DENSITY 

In the past several decades, QCD sum rule has been widely used to study the hadron structures and proven to be 
a very powerful non-perturbative method [47l - l49j . We consider the two-point correlation function: 

n^-vCg') - i [ d^xe''J-{0\T[J^{x)Jt{0)]\0) 



= -ni(,2)(.9,.-Mf:) + no(g^)^, (5) 

There are two independent parts of II^j^ with different Lorentz structures because is not a conserved current. 
Ili{q^) is related to the vector meson while Ilo{q^) is the scalar current polarization function. 

At the hadron level, the correlation function is expressed by the dispersion relation with a spectral function: 

ni(g^)= r . , (6) 

where the lower limit of integration is the square of the sum of the mass of all current quarks (omitting the light 
quark mass). The spectral function is defined as: 

pis) = Y,S{s^rnl){0\r^\n){n\rj+\Q) 

n 

= /^i5(s — m^) + continuum, (7) 
where the usual pole plus continuum paramctrization of the hadronic spectral density is adopted. 
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On the other hand, the correlation function can also be calculated at the quark-gluon level via the operator product 
expansion(OPE) method. Using the same technique as in Refs. [13, SH Ey, [s^], we evaluate the Wilson coefficient 
up to dimension eight while the coordinate space expression for the light quark propagator and the momentum space 
expression for the charm quark propagator are adopted: 

•ra& • \n -1 Slab Xab 2 

where x = ■y^ix'^, p = 7^p^, (qgsCrGq) = {gsqcr^^Gfj^^q) , {g^GG) = {glG^^G^^'^), a and b arc color indices. The 
dimensional regularization is used throughout our calculation. The spectral density is then obtained with: p{s) = 
ilmn(g2). 

One of the important assumption in the QCD sum rule approach is the quark-hadron duality, which ensures the 
equivalence of the correlation functions obtained at the hadron level and the quark-gluon level. To improve the 
convergence of the OPE series, the Borel transformation is performed to the correlation functions at both levels. 
Considering the spectral function defined in Eq. ([T]), we arrive at: 

/•So 

f2^^-rnl/M% ^ / dse-^/^^« , (9) 

J 4m2 

where sq is the threshold parameter. Then we can extract the mass mx- 



(10) 



The spectral extracted densities of all the tctraquark currents in Eq. (|3]) are listed in the Appendix. For each 
quantum number, we just list the expressions of four currents. Others could be obtained conveniently by simple 
replacement rric — > —rric. We collect the terms that are proportional to the light quark mass ruq into the expressions 
of the spectral densities. These terms give small contributions to the correlation functions involving the qcqc system 
and hence can be ignored. In the case of the scsc system, however, they give important corrections because the strange 
quark mass is much larger than m„ and m^. Both the quark condensate {qq) and quark gluon mixed condensate 
{qgsCrGq) appear with both and nic, which is very different from the case in the pseudoscalar channel [4^, where 
only the nig related terms contribute. We neglect the three gluon condensate g^{fGGG) because they are strongly 
suppressed and negligible [46j . 



IV. QCD SUM RULE ANALYSIS 

We perform the QCD sum rule analysis using the following parameter values of the quark masses and various 
condensates [irl lsTI - isi l : 

m„(2 GeV) = (2.9 ± 0.6) MeV , 
md(2 GeV) = (5.2 ± 0.9) MeV , 
771,(2 GeV) = (4.0 ± 0.7) MeV , 
777,(2 GeV) = (lOlla?) MeV, 
777c(777c) = (1.23 ±0.09) GcV, 
7776(7776) = (4.20 ±0.07) GeV, 
(qq) = -(0.23 ±0.03)3 GeY^ , 
(qgsCrGq) = -MQ{qq) , 
Ml = (0.8 ± 0.2) GeV^ , 

(ss)/(gg) = 0.8±0.2, (11) 
{glGG) = 0.88 GeV^ . 

where the up, down and strange quark masses are the current quark masses in a mass-independent subtraction scheme 
such as MS at a scale fi = 2 GcV. The charm and bottom quark masses are the running masses in the M S scheme. 
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A. The tetraquark systems with J — I and 1 



The stability of QCD sum rule requires a suitable working region of the threshold value sq and the Borel mass AIb- 
In our analysis, we choose the value of Sq around which the variation of the extracted mass mx with M"^ is minimum. 
The working region of the Borel mass is determined by the convergence of the OPE series and the pole contribution. 
The requirement of the convergence of the OPE series leads to the lower bound M^j„ of the Borel parameter while 
the constraint of the pole contribution yields the upper bound of Mg. 

We first study the interpolating currents with J^^ = 1 and 1 in the qcqc systems. For these currents, the 
contribution of the four quark condensate {qq)'^ is negative and its absolute value is bigger than other condensates in 
the region of Mg < 3.1 GeV^. It is the dominant power contribution to the correlation function in this region. In 
fact, the quark condensate (qq) in the spectral density of the 1 currents is proportional to the light quark mass 
niq and vanishes if we take rriq = 0. For the 1 ^ currents, we require that the perturbative term be larger than five 
times of the four quark condensate to ensure the convergence of the OPE series, which results in the lower bound of 
the Borel parameter, M^j„ ^ 2.9 GeV^. The constraint for the currents is stricter due to the vanishing (qq) 
condensate. The pole contribution (PC) is defined as: 



PC 



(12) 



It is clear that the PC depends on the threshold value sq. As mentioned above, sq is chosen to ensure the minimum 
variation of mx with A/^. For example, we choose so ~ 26 GeV^ for Jg^ with J^'^ = 1 ^ in Fig. [l}i. By requiring 
the pole contribution be larger than 40%, we obtain the upper limit M^^^ of the Borel parameter. 
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FIG. 1: The variation of mx with so (a) and M%{h) corresponding to the current Je^ for the 1 qcqc system. 
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FIG. 2: The variation of rrfx with so (a) and Mg(h) corresponding to the current Je/j for the 1 ^ scsc system. 
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After performing the QCD sum rule analysis in the working range of the parameters obtained above, only the 
currents Jg^^, Jt^ and Jg^ with J^'^ = 1 have the stable mass sum rules. In Fig. [TJ we show the variation of mx 
with the threshold value sq ^nd Borcl parameter Mg for the current Jq^ in qcqc system. The plateau in the region of 
10 ^ 14 GeV^ in Fig. [T^ is just an unphysical artifact because both the numerator and denominator in Eq. ([TU]) are 
negative within this region. The situation also occurs in the pseudoscalar channel in Ref. [46| . The variation of mx 
with the Borel parameter is weak around the region sq ~ 5.1^ GeV^, as shown in Fig. [TJd. For Ji^j, J2^j, Js^^, JA^^ 
and J^fj,^ the stability is so bad that the extracted mass mx grows monotonically with the threshold value sq and the 
Borcl parameter Mb- These currents may couple to the 1 ^ states very weakly and the continuum contribution may 
be quite large, leading to the above unstable mass sum rules. 

We show the Borel window, the threshold value, the extracted mass and the pole contribution corresponding to the 
tetraquark currents Jq^ ^ Jg^ with J^*^ = 1 ^ in the qcqc system in Table HIl The results of the 1 system are 
listed in Table Hill We present the numerical results for the currents which lead to the stable mass sum rules in the 
working range of the Borel parameter. Only the errors from the uncertainty of the threshold values and variation of 
the Borel parameter are taken into account. Other possible error sources include the truncation of the OPE series 
and the uncertainty of the quark masses, condensate values and so on. 





Currents so(GcV2) 


[A^^nin'^^maxKGeV^) mx(GeV) 


PC(%) 






5.1^ 


2.9 ~ 


' 3.9 


4.67 ±0.10 


50.2 


qcqc system 




5.22 


2.9 - 


' 4.2 


4.77 ±0.10 


47.4 






4.92 


2.9 - 


' 3.4 


4.53 ±0.10 


46.3 




Jin 


5.02 


2.9 - 


' 3.4 


4.67 ±0.10 


44.3 






5.02 


2.9 - 


' 3.4 


4.65 ±0.09 


45.6 




Jin 


4.92 


2.9 - 


' 3.3 


4.54 ±0.10 


44.4 






5.12 


2.9 - 


' 3.7 


4.72 ±0.09 


44.8 


scsc system 


JbiJ. 


5.02 


2.9 - 


' 3.6 


4.62 ±0.10 


42.8 




Jen 


5.32 


2.9 - 


' 4.3 


4.84 ±0.10 


47.3 






5.32 


2.9 - 


' 4.3 


4.87 ±0.10 


46.2 




Jifi 


5.22 


2.9 - 


' 4.1 


4.77 ±0.10 


44.1 






11.02 


7.2 ~ 


' 8.6 


10.53 ±0.11 


44.2 


qbqb system 




11.02 


7.2 - 


' 8.6 


10.53 ±0.10 


44.1 




J&n 


11.02 


7.2 - 


' 8.6 


10.49 ±0.11 


44.7 




Jiti 


11.02 


7.2 ~ 


' 8.1 


10.62 ±0.10 


41.2 






11.02 


7.2 ~ 


' 8.4 


10.56 ±0.10 


43.8 


qbqb system 




11.02 


7.2 - 


' 8.3 


10.63 ±0.10 


42.4 






11.02 


7.2 ~ 


' 8.3 


10.62 ±0.09 


42.5 






11.02 


7.2 ~ 


' 8.3 


10.59 ±0.10 


43.1 



TABLE II: The threshold value, Borel window, mass and pole contribution corresponding to the currents with J^'^ = 1 in 
the qcqc, scsc, qbqb and sbsb systems. 

The analysis can easily be extended to the scsc systems. We keep the mg related terms in the spectral densities. 
These terms give important corrections to the OPE series and are propitious to enhance the stability of the sum rule. 
Especially for the currents with J^'~' = 1 . (ss) is now proportional to the strange quark and larger than (ss)2 
in the Borcl window, which is very different from the qcqc system. Using the parameters in Eq. (jl2p . we also collect 
the numerical results for the 1 and 1 scsc systems in Table HIl and Table |lTT] respectively. Obviously, the scsc 
systems have better stabilities than the qcqc systems. We show the variation of mx= with sq and for the current 
Jefj, with J^"^ = 1"+ in Fig. H which is very similar to Fig. [Tj except the value of sq around which the variation of 
mx with Mb is minimum. The extracted mass of the scsc state is a little higher than that of the qcqc state. The 
mass difference mx= — mx for the same interpolating current is about 0.2 GeV, which is about 2(tos — m,) within 
the errors. 

The extracted mass of the qcqc state with J^'~'' = 1^^ in Table Hill is about 4.6 ^ 4.7 GcV, which is consistent with 
the mass of the meson y(4660). One may wonder whether y(4660) could be a tetraquark state. 

Properties of the bottomonium-likc analogues tend to be very similar because of the heavy quark symmetry. Re- 
placing mc with mt in the correlation functions and repeating the same analysis procedures done above, we collect 
the relevant results of the qbqb and sbsb systems in Table [ll] and Table IIIII One notes that the bottomonium-like 
systems are less stable than the corresponding charmonium-like systems. 
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Cm-rents 


.so(GeV") 




mx(GeV) 


PC(%) 






5.0^ 


2.9 ~ 


3.6 


4.64 ±0.09 


44.1 


qcqc system 




5.02 


2.9 - 


3.6 


4.61 ±0.10 


46.4 






5.22 


2.9 - 


4.1 


4.74 ±0.10 


47.3 






0.4 




4.0 


/I 09 -1- n 1 n 


OU.O 






5.02 


2.8 - 


3.5 


4.64 ±0.09 


48.6 


scsc system 




4.92 


2.8 - 


3.4 


4.52 ±0.10 


45.6 






5.42 


2.8 - 


4.5 


4.88 ±0.10 


51.7 






5.3 


2.8 - 


4.3 


4.86 ± 0.10 


46.0 






4.8 


2.8 - 


3.1 


4.48 ± 0.10 


43.2 


qbqb system 




11.02 


7.2 - 


8.5 


10.51 ±0.10 


45.8 




7, 


1 1 n2 
11. u 


7 
1 .z ~ 


C Q 
0.0 


lU.DU It U. lU 


A7 n 
4 i .u 




" z/i 


11.02 


7.2 ~ 


8.4 


10.55 ± 0.11 


43.6 


sbsb system 




11.02 


7.2 - 


8.4 


10.55 ±0.10 


43.7 




Jifi 


11.02 


7.2 - 


8.4 


10.53 ±0.11 


44.3 






11.02 


7.2 - 


8.2 


10.62 ±0.10 


42.0 






11.02 


7.2 - 


8.4 


10.53 ±0.10 


44.1 



TABLE III: The threshold value, Borel window, mass and pole contribution corresponding to the currents with J = 1 in 
the qcqc, scsc, qbqb and sbsb systems. 



B. The tetraquark systems with J 



and r 



We study the cmrents with J^'~'' = 1++ and 1^ in this subsection. The spectral densities of these currents are 
very similar to that of the 1 ^ and 1 currents, as shown in the Appendix. The analysis shows that the OPE 
convergence becomes worse than that in the vector channel. In this channel, the quark condensate {qq) is larger than 
any other condensates for all the currents. The OPE convergence of the currents J^^, J^f^, Jjf^, Js/i is a little better 
than that of Ji^, Jj^, Jg^, J4^. 

Only the currents Ja^ and J^^^ in the 1"^"*" qcqc system display stable mass sum rules. We obtain the working region 
of the Borel parameter in 3.0 < < 3.4 GeV2 while taking sq ~ GcV^ for current J4^. The variations of mx 
with the threshold value Sq and Borel parameter are shown in Fig. [31 from which the Mg dependence is very 
weak around the chosen threshold values. Taking into account only the errors from the variation of Mb and sq, the 
extracted mass is mx = 4.03 GeV, which is slightly above the mass of X(3872) within the errors. 

The scsc, qbqb and sbsb systems can be studied conveniently by replacement of the parameters, including the quark 
masses and the various condensates. The numerical results are listed in Table Hvl for the 1++ systems and Table IVl 
for the l"* systems. Now the bottomonium-like systems are more stable than the corresponding charmonium-like 
systems. 
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FIG. 3: The variation of mx with so (a) and M^(b) corresponding to the current Jn^ for the 1^"*" qcqc system. 
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Currents 


so(GeV") 




mx(GeV) 


PC(%) 


qcqc system 




4.6^ 


3.0 - 


' 3.4 


4.19 ±0.10 


47.3 






4.5 


3.0 - 


' 3.3 


4.03 ± 0.11 


46.8 


scsc system 




4.62 


3.0 - 


' 3.4 


4.22 ±0.10 


45.7 






4.52 


3.0 - 


' 3.3 


4.07 ±0.10 


44.4 




J^ti 


10.92 


8.5 - 


' 9.5 


10.32 ±0.09 


47.0 


qbqb system 


Jin 


10.82 


8.5 


' 9.2 


10.22 ±0.11 


44.6 




J7n 


10.72 


7.8 - 


' 8.4 


10.14 ±0.10 


44.8 




Jsn 


10.72 


7.8 - 


' 8.4 


10.14 ±0.09 


44.8 




J^ij. 


10.92 


8.5 ~ 


' 9.5 


10.34 ±0.09 


46.1 


sbsb system 




10.32 


8.5 - 


' 9.1 


10.25 ±0.10 


43.3 




J Til 


10.82 


7.5 - 


' 8.6 


10.24±0.11 


47.1 






10.82 


7.5 ~ 


' 8.6 


10.24 ±0.10 


47.1 



TABLE IV: The threshold value, Borel window, mass and pole contribution corresponding to the currents with J = 1^"*" in 
the qcqc, scsc, qbqb and sbsb systems. 





Cm-rents so(GeV^) 


[A^^i,,^^f^ax](GeV^) 


mx(GeV) 


PC(%) 






4.62 


3.0 - 


' 3.4 


4.16 ±0.10 


46.2 


qcqc system 


Jin 


4.52 


3.0 - 


' 3.3 


4.02 ±0.09 


44.6 




Jbn 


4.52 


3.0 - 


' 3.4 


4.00 ±0.11 


46.0 




Jen 


4.62 


3.0 - 


' 3.4 


4.14 ±0.09 


47.0 




Jsfi 


4.72 


3.0 - 


' 3.6 


4.24 ±0.10 


49.6 


scsc system 


Jin 


4.62 


3.0 - 


' 3.5 


4.12 ±0.11 


47.3 




Jbn 


4.52 


3.0 - 


' 3.3 


4.03 ±0.11 


44.2 




Jen 


4.62 


3.0 - 


' 3.4 


4.16 ±0.11 


46.0 




Jsn 


10.62 


7.5 ~ 


' 8.5 


10.08 ±0.10 


45.9 


qbqb system 


Jin 


10.62 


7.5 ~ 


' 8.5 


10.07 ±0.10 


46.2 




Jbn 


10.62 


7.5 ~ 


' 8.4 


10.05 ±0.10 


45.3 




Jen 


10.72 


7.5 ~ 


' 8.7 


10.15±0.10 


47.6 




Jsn 


10.62 


7.5 - 


' 8.3 


10.11 ±0.10 


43.8 


sbsb system 


Jin 


10.62 


7.5 ~ 


' 8.4 


10.10±0.10 


44.1 




Jbn 


10.62 


7.5 ~ 


' 8.3 


10.08 ±0.10 


43.7 




Jen 


10.72 


7.5 ~ 


' 8.5 


10.18±0.10 


46.5 



TABLE V: The threshold value, Borel window, mass and pole contribution corresponding to the currents with J^*^ — l'^ in 
the qcqc, scsc, qbqb and sbsb systems. 

V. DECAY PATTERNS OF THE CHARMONIUM-LIKE STATES 

The decay properties of the charmonium-like states are important for the study of their structures and detections 
at experiments. In this section, we study the decay patterns of the charmonium-like states with J^'^' = 1 '", 1 , 1++ 
and 1"' . Only the two-body hadronic decay is considered. Replacing the D meson by the B meson, one gets the 
decay patterns of the bottomonium-like states so long as the kinematics allows. 

The G-parity of a charmonium-like state is defined as G = C ■ (—1)''^, where / is the isospin. By considering the 
conservation of the angular momentum, P-parity, C-parity, isospin and G-parity, we collect the S-wave and P-wave 
decay modes of these charmonium-like states in Tabic IVII and Table I VIII For the vector channel, one notes that the 
S-wave decay modes arc dominant and the final states always contain a S-wavc and P-wavc meson pair. Such a decay 
pattern is also speculated to be characteristic of the hybrid meson. In fact, the tctraquark state mixes easily with 
the hybrid state cGc. In quantum field theory the charmonium-like tetraquark operator and the cGc hybrid operator 
with the same quantum numbers probably couple to the same physical state. In Table IVlTl the S-wave decay modes 
J/ifjui and J/ipp are hsted in the l"*""*" channel, which is consistent with the decay properties of X(3872) [l,!!^. Up 
to now, no experimental signals are observed for the charmonium-like 1 and 1"' states. y(4360) and F(4660) are 
only observed in the ip{2S)'K~^Tr~ channel [H, |5l|. The possible decay modes listed in Table IVll and Table IVlTl may 
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be useful to the future search of these interesting charmonium-hke and bottomonium-like states at the experimental 
facilities such as Super-B factories, PANDE, LHC and RHIC. 

ju jfu S'-wave P-wave 

(2007)"Z3(5 (2400)^ + C.C., L»^'(1865)i:'''(1865), i:)*(2007)"i:»*(2007)^ 

0- 1— A(2420)°50(i865) + C.C., Xco{lP)hi[inQ) , Xci{lP)hi{inQ) , 

Z?i(2420)05*(2007)0 + c.c.,J/V/o(980), J/^r/, i>{2S)i , ^(25)?7, ric{lS)uj,r^,{2S)uj 
V'(25)/o(980), Xco(lP)^, Xci(lPV, 

£)*(2007)°Z)S(2400)° + c.c, D°(1865)i5°(1865), D*(2007)°L'*(2007)°, 

1+1— i:)i(2420)''L>"(1865) + c.c, Xco(l-P)&i(1235), Xci(l-P)6i(1235), J/ipn, 

i^i(2420)"5* (2007)0 + c.c, Xco(l-P)p, ^'(2S)7r, ?7,(lS)p, 7,,(25)p 

Xci(lP)p, ai(1260)J/^-, 6i(1235)77e(15), 

L)*(2007)'JZ3S(2400)'J + CC, i:»"(1865)i:'"(1865), i:)*(2007)"i:»*(2007)^ 

0+1-+ £»i(2420)"50(1865) + CC, Xco(lP)/o(600), Xco(lP)/o(980), 

£)i(2420)OL»*(2007)° + c.c, rjc{lS)T],r]AlS)r]' , J/^juj, il;{2S)w 

/i(1285)ry,(15), Xci{lP)v, Xci{lPW Xci(lP)/o(600), Xci(lJ^)/o(980), 

L>*(2007)°I)*_(2400)° + CC, ai (1260)77^(15'), i:»°(1865)i!»°(1865), D* (2007)0/)* (2007)°, 

1- 1-+ Di(2420)OL)0(1865) + c.c, 6i(1235) J/tA, Xco(lP)ao(980), r?c(15)7r, Vc{2S)tt, 
£>! (2420)°!)* (2007)" + c.c, Xci(l-P)7r J/^p, ^{2S)p 

TABLE VI: The possible decay modes of the 1 and 1"+ charmonium-like states. 





J^^' S'-wave 


P-wave 


0+ 


1++ i:»"(1865)i:'*(2007)" + cc, Xci(l-P)/o(600), 


Xco{iPH XciiiPH 




i:'(1870)+i:>*(2010)- -hcc, J/i)i^j 


?7,(lS)/o(600) 


1- 


1++ i:>O(1865)i5*(2007)° -Kcc, 


Xco(l^)^, 




L'(1870)+i:»*(2010)- 4- cc, J/ijjp 

1-1— nU^occN n*^ooo-7\0 , „„ 


Xci{lP)n, ?7,(lS)ai(1260) 

„ M o\u /'II -TAN T/^i.-C rc:r\rt\ 



D{1870)+ D* (2010)- + C.C, Xcoi^Pp 
J/ip-q, J/i'i, ric{lS)uj 

1+1+- D0(1865)L'*(2007)° + CC, J/V'tt, iI:{2S)t:, 7/c(1S')6i(1235), Xco(l^')p 
i:>(1870)+i:>*(2010)- + c.c, T]c{lS)p 

TABLE VII: The possible decay modes of the 1++ and 1"' charmonium-like states. 



VI. SUMMARY 



We have constructed a matrix O which is composed of the tetraquark operators with different Lorentz structures. 
The charge-conjugation transformation of the matrix is equal to its transpose. With this relation, we can define 
the symmetric matrix S and antisymmetric matrix A. Considering the color structure, the elements of S and A are 
the tetraquark operators with definite C-parities. Then we can obtain all the tetraquark interpolating currents with 
J^C" = 1-+, 1 — , 1++ and 1+-, as shown in Eqs. ©-Q- 

At the hadronic level, there exists big difference between the molecular states and tetraquark states. The molecular 
states are commonly assumed to be bound states of two hadrons formed by the exchange of the color-singlet light 
mesons while the tetraquark states are generally bound by the QCD force at the quark gluon level. However, within 
the framework of the QCD sum rule approach, the only difference lies in the interpolating current used in the study of 
the molecular and tetraquark states. Other procedures including the operator product expansion, the calculation of 
the Wilsion coefficient and numerical analysis are the same. In principle, if we exhaust all the possible molecular-type 
currents and all the possible tetraquark-type currents, we can rigorously show that these two sets of interpolating 
currents are equivalent by using the Fierz rearrangement. 

However, there exists important difference between one single molecular- type current and one single tetraquark-type 
current. For example, every single tetraquark-type current is a linear combination of several independent molecular- 
type currents. In this respect, one well-known example is the light scalar-isoscalar sigma meson. The tetraquark-type 
current (or their combination/mixing) leads to a better result than the simple pion-pion molecular current. It's 
possible to distinguish the tetraquark and molecular structures after exhaustive and comprehensive hard work, which 
is one of the motivations of our present systematical investigation. 
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By studying the two-point correlation functions, we have calculated the spectral densities of these currents at the 
quark-gluon level. The four quark condensate {qq)"^ is dominant for all the currents with J^'-' = 1 and 1 in 
the qcqc systems. For the currents with J^'~' = 1++ and 1^ , however, the quark condensate (qq) becomes the most 
important corrections to the correlation functions. These properties of the spectral densities lead to a better OPE 
convergence for the currents in the vector channel than that in the axial-vector channel. The related terms are 
kept in the spectral densities in order to study the contribution of the strange quark in the scs(5 system. Actually, 
they give important corrections and enhance the stabilities of QCD sum rules. 

The tetraquark assignments of X{3872) and F (4660) have been studied within the framework of the QCD sum rule 
approach, as mentioned in the introduction. In Ref. [17[, the current was used to study the 1+"'" qcqc system with 
the extracted mass around 3.92 GeV. The current Ja^ was used to study the 1~~ scsc system with the extracted 
mass around 4.65 GeV Ref. [l^. One difference of our present analysis and the previous ones lies in the criteria of 
fixing the Borel window and the value of the threshold parameter sq, which leads to the slightly different extracted 
masses of the states. We have imposed a strict requirement that (1) the pole contribution be larger than 40% and 
(2)dual stability, i.e., the variation of the extracted mass with both so and the Borel parameter be minimum. The 
other difference is that we have exhausted the tetraquark interpolating currents. 

In the working range of the Borel parameter, only the currents Ji^^Ji^ and Jj^ with J^'-^ = 1 display stable 
QCD sum rules in the qcqc system. The extracted mass is around 4.6 ~ 4.7 GcV from these currents, which is 
consistent with the mass of the meson F(4660). This result implies a possible tetraquark interpretation for y(4660). 
In the scsc system, all currents except J5^, Jg^ have stable QCD sum rules. The mass difference tox" — fnx ~ 0.2 
GeV for the same type of the interpolating current, which is roughly 2{rns — mq). For the 1 bottomonium-like 
states, the masses lie around 10.5 GeV and 10.6 GeV for the qhcjb and shsi systems, respectively. The Borel window 
for the currents in the axial- vector channel is very small because of the bad OPE convergence. For the currents with 
J^*-^ = l"*""*" in the qcqc system, only and lead to stable QCD sum rules. The same situation occurs in the scsc 
system. The extracted masses are about 4.0 ^ 4.2 GeV, which is 0.1 ^ 0.2 GeV higher than the mass of X(3872). 
The masses of the 1++ bottomonium-like states are about 10.2 GeV for both the qh(jb and shsh systems. 
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Appendix A: THE SPECTRAL DENSITIES 



In this appendix wc show the spectral densities of the tetraquark interpolating currents defined in Eq. Various 
power corrections include the four quark condensate (qq)'^, quark gluon mixed condensate (qgsCrGq) and dimension 
eight condensate {qgsCrGq){qq): 

p{s) ^ pP'^'-^s) + p<«'?>(5) + p(GG)(^) _^ p(9'?)'(s) + (^^^ 

The integration limits in the expressions are: 

l + y/l- Aml/s _l- ^Jl- Aml/s 

« - 1 _ R - 

Pniax ~ J- Pmin — o • 

as — 
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1. The spectral densities for the currents with J^*^ — 1 
For the interpolating current Ji^: 



Pi {s) 



Pi 



{qgsCrGq) 



is) 



p'r'\s) 



3 X 287r6 



M 

87r4 



da I dp 

Ip. 



da j dp 



{l-a-P)[{a + l3)ml-af3s] ^ / 3(1 + a + l3)[{a + p)ml " "/^s] 
f '"^-(^"""^^' +24...n.,(l-a-/3)L 



/3 



[(a + ^)m^ - a/3s] / 2mc(l - a - /3)[(a + /S)™^ _ q,^<, 
32 X 2i27r6 



J(a + /3 - 3)m2 - 4a/3s] ^ , 



da 



dp 



{I -a- Pfmf 



[96/32 _^ 3g^^ _ (1 „ a - /3)(5a + 48/3)][(a + P)ml - a/3s] 30[(a + P)ml - a/3s] 



a/32 



/32 



96[2(a + /3)m2 - 3a;3s] + 16(1 - a - /3)m2 



a 



48(1 - a - /3)[(a + P)ml - a/3s]s 



a 



3 X 277r4 



I da 1 


°'d/3| 







TOj48a/3 + (1 - a - /3)(7a + 6/3)] [(a + /3)m2 - a/3s] 
2mc(l - a - /3)[2(a + /3)to2 „ 3^^/35] Qmq[{a + /3 - 2)m2 - a/3£ 



/?2 



{qq)'^{8ml ~ 6mcmq + s) / ^ , 

36^2 4m,/ i 



n 



'(Af 



127r2 



2m'; m2 A/rQ; 



(A2) 



For the interpolating current J2^: 



p'rks) 



(qgscGq) 
Pi 



pr (.), pt\s) = (s), (^) = :^pr (^) , 



(9-9)^ 



3 {qq? 



32 X 2i27r6 



(1 - a - Pfml 



da / d/3 



(1 - a - ;3)(13a + 72;3)[(a + ^)m2 - a/3s] 42[(a + ;3)m2 - a/3s 



24[(1 + 5a + 5/3)m2 - 12a/3s] 



a 
2 



a/32 

12[(a + /3)m2 - a/35 



■[(5 - a - /3)m2 - (1 - a - pfs - 2a/3s] ^ , 



mq{qgsaGq){8m- + s) /- ^ , 

64^^ v/l-4m2As 

r 3n „ ^ _ ^usj-^ _,_ ^)^2 _ 5^,^^! 20[(a + P)ml ~ a^s 
dpi 



.{qgsCrGq) 



27^4 



da 



13' 



n 



(qgs<yGq)(qq) 



{qgs(jGq){qq) 
87r2 



da{ A/|(l - a) 



1 - a 



a2A/2 



(A3) 
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For the interpolating current J^^: 



(1 - a - /3)2m2 
8(l-a-/3)m2 



da 



/3„ 



dp 



[36a/3 + 24^(1 - a - P) - a{5 + a + /3)][(a + /3)m2 - a/3s] 48[(a + ^)m2 - 2a/3s 



a/32 

48(1 - a - /3)[(a + /3)m2 - a/3s]s 



a 



1927^4 ^""^/^ 



3 X 277r4 



da 



d/3 



6mc(l — 5q! — /3)[(q! + 13)^1 — aj3s 



TOc(l - a - I3)[l'b{a + /3)to2 - 29a^s] 6mq[(a + /3 - 2)to2 - a/J, 



/32 



+ 



n 



{qgs<yGq){qq) 



{Ml) 



{qgsCFGq){qq) , /.„ o n , r2 
48^^ da|(2-3a)A/^ 



2m? 



4m4 
a2A/2 



For the interpolating current 



4ti- 



(A4) 



pf'-*(.) = 3pr\s), pr is) = 3pr (^), (^) = (^) , 



peri / 



Aqq) 



(99) / 



Aqqf 



32 X 2i27r6 



r-/3™o. r (1 _ Q, _ ^)2^2 r48[m2(l + 5a + 5/3) - 12a/3s] 
da / d/3< 



a 



+ 



[216a;3 + (1 - a - /3)(65a + 144;3)][(a + /3)m2 - a/3s] 210[(a + (3)ml - a/3s 



a^2 



/32 



12[(a + /3)m2 -a/3s] 



[5(5 - a - ^)m^ - 5(1 - a - l3)^s ~ 24(1 - a - /3)/3.s - 10a/3s] 



mg{qgsaGq){8ml + .s) 
32^i^ 



^/l - 4m2/i 



2^714 



/-/3™ax rmr[28a^-(23a+12/^)(l-a-/^)][(a + /3)m2-a;3s] 

"1.. 

22toc(1 - a - /^)[2(a + /3)to2 - 3af3s] I2mg[{a + /3 - 2)m2 - a/3s] 
^2 



n 



(qgsiyGq)(qq) 



{Ml) 



{qgsaGq){qq ) 
87r2 



da<^ A/^(2-3a) - 2m^ 



a2i\/2 



(A5) 



From these results the expressions for the currents Js^, Jq^,, Jj^ and Js^ can then be obtained conveniently by the 
replacement rric — !> —mc- 



Pi{s) 



rnc^ — trLc 



> P5{s),P2{s) '"""^ P8{s),P3{s) 



rfic—^ — mc 



> P7(s),P4(s) 



> P6(s) . 



(A6) 
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2. The spectral densities for the currents with J^'~^ — 1 
For the interpolating current Ji^: 



(GG), X 

Pi (s) 



Pi 



p'^'Us) 



1 



3 X 287r6 



da 



dp 



{l-a~P)[{a + p)ml~aps] 



3(1 + a + I3)[{a + (3)ml - afis] - 2ml{l - a - /3)H , 



I. 

{glGG) 

32 X 2i27r6 



da 



/3„ 



d/3 



[(a + I3)ml - a/3s][(5 - a - I3)ml + 2a/3s 



/3„ 

da 



dp 



48(1 -a 5(5 + a + /?) 



36 



(1 - a - pf[{a + /3)mg - a/gsjm^ 
a2/3 

16(1 - a - /3)2[(1 - 7a - ll3)ml + 12aPs]ml 48(1 - a - /3)[(a + ^)m2 - a;3s]s 



(qg.aGq), . _ mq{qg,,<7Gq) {16ml - s) 



3 X 257r4 
(1 - a - P)mc 



3 X 277r4 



da 



/3 

(g(?)2(16m2 - s) 
367r2 



6[{a + I3)ml - al3s] 27(a + /3)m2 - 49a/3s 
a ^ 



dp 

6m,((2 + a + /3)m2 - a/3s) 



a 



^1 - 4m2/5 



2(1 - 2a)?7i2 8m4 



a^Afl 



For the interpolating current 

pert 

{glGG) 



pr\s) 



(1 - a - ^)2[(a + P)ml - aps] [72(1 - a - /3)m2 (29 + 13a + np)ml ~ UaPs 



a^P 



a 



P 



(A7) 



12[(a + P)ml - aPs][{5 -a- P)ral + 2a/3s] 24(1 - a - /3)2[(1 - 7a - 7P)ml + 12aPs]ml 



aP 



mg(qgsCrGq){16ml - s) 



647r4 



/- , . niciqgsCrGq) 



d/3 



(l-a-/3) 



12[(a + /3)TO2 _ Q,^s] 3(a + /3)m2 - 13a/35 



a 



/3 



n 



{qgaCrGq){qq) 



(M, 



{qg^aGq) {qq) 
87r2 



(1 - a)M| + 



(1 — a)ml 2m^ 



a2M2 



(A8) 
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For the interpolating current J^^: 



{GG), X 



{glGG) 



da 



dp 



32 X 2i27r6 

(1 - a - /3)2[(a + P)ml - aPs]ml [24(1 - a - (3) 5 + a + P 



a 



36 



8(1 - a - /3)2[(1 - 7a - 7/3)m2 + 12a/3s]m2 48(1 - a - /3)[(a + /3)to2 - aPs]s 



■mq{qgsuGq){lQml 
3 X 267r4 

(1 - a - ^)mc 



-yl — 4mf7s 



3 X 277r4 



da 



6[(a + /3)m2 - a/3s] 9(a + /3)m2 - 19a/3s' 
a ^ ^ 



d/3 

6TOg((2 + a + I3)ml ~ a/3s) ' 



n 



(qgs<yGq){qq) 



{Ml) 



{qgsCrGq){qq) 



487r2 



da<^ (2 - 3a)A/| 



2(2 - a)ml 4m^ 



For the interpolating current J4^: 



p(jl9.<rGq)^^^^ 



da 



'iffr\sl pi'^'Hs) ^ 3pf'^>(.), > (s) = 3p<'"'> is) 

{glGG) 

32 X 2i27r6 

(1 -a-/3)2[(a + /3)m2 -a^s]m2 ri44(l-a-/3) 65(1 - a - /3) 210 



dp 



I3„ 



a^(3 



216 



60(5 - a - ^)[(a + (3)ml - a/3s]m2 48(1 - a - /3)2[(1 - 7a - 7/3)m2 + 12a/3s]m2 

fl ' 3 

ap a'' 

288(1 - a -/3) 60(l-a-/3)2 



a ap 
mq{qgsCrGq){16ml - s) 



120 



[{a + P)m1 — aPs\s 



327r4 

(1 - a - P)mc 



a/1 - 4m2/s + 



{qgsCrGq) f" 



da 



3 X 2^774 

24[(a + /3)m2 - a/3s] 87(a + P)ml - 161a/3s' 



dp 



(A9) 



36m5((2 + a + /3)to2 _ 0,^35) ' 



a 



j,<^9.^0qH'^q) ^Ml) = (^Ss^^f'il) f'dal{2-3a)Ml + ^-^^^^-4^}e °<-°)-^ . (AlO) 

8^ Jo i a^Mg 



From these results the expressions for the currents J^^, Jg^, J^^^ and Jg^ can then be obtained conveniently by the 
replacement nic —ttIc'. 

Pl[S) > P5{s),p2(s) > ps{s),P3{s) > p7(s),P4(s) > P6(S) . (All) 
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3. The spectral densities for the currents with J^^ = 1"*""*" 
For the interpolating current Ji^: 



Pi [s] 



da 



dl3- 



+ 12mcm,(l -a- l3)[{a + l3)ml - 3a/3s] [ , 



{1 + a + I3)[{a + I3)ml ^ al3s] 



M r 

87r4 



da 



dp 



[{a + I3)ml - a/Ss] f TOc(1 -a- /3)[3(a + I3)ml - la/Ss] 
ap' 

32 X 2i27r6 



[(4 + a + I3)ml - 3a/3s] 



da 



dp 



{l-a~ P)^m^ 



96[{a + P)ml - 2aPs] 5(5 + a + ip)[{a + P)ml - a^s] 



a 



&[{a + P)ml -aps] 



aP 
87r4 



4(l-a-/3)(5m2-2a/3s) 



P' 

5ml{l + a + P) 



a~ P)[i{a + P)ml- aPs] I2mqml {29171^ + 5mq)[3{a + P)ml - 5aPs] 



P' 



P 



P 



mc{qq) 



12^2 (4me + 2m, - - Aml/s , 



n 



(qgs<yGq)(qq) , , r2 



{hi 



Bl 



32 X 257r2 



i{l - a) a^Ml 



(A12) 



For the interpolating current J; 



2tJ,- 



32 X 2i27r6 



(1 - a - Pfml 



a'' 



48[(a + /3)m2 - 2a/3s] (5 + a + 7/3)[(a + P)'ml ~ aps\ 



6[(a + /3)m2 - a^s] 



aP 

mlmq{qgsaGq) 
167r4 



4(l-a-/3)(m2 + 2a/3s) 



/32 

-m2(l + a + /3) 



,(l-a-^)[15(a + /3)m2 -29a/3s] Umgml {Urric + mq)[3{a + P)ml - 5aPi 



/32 



/3 



/3 



n 



{qgs<jGq){qq) 



(Ml) 



{qgsCrGq){qq) 
32 X 2^712 



V r„ , . 2(6-5a)TO2 

dai 3aMl + -^ '—^ 

I Q;(l-a) 



24m4 



(A13) 
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For the interpolating current Js^: 



(GG) , ^ 



1 



3 X 287r6 



da 



dp 



(1 - a - /3)[(a + /3)to2 _ a^sf / 9(1 + a + + I3)ml - ajSs 



4m^(l - a - /3)(5 + a + + /3)ml - ajBs 



+ 36TOcmg(l — a — /3)[(a + /3)m^ — 3q!/3s] > , 



M r 

87r4 



da 



dp 



[{a + j3)ml- a/3s] / 7/1^(1 - a - /3) [3(a + /3)to^ _ Jq,/? . 
32 X 2i27r6 



+ mj(14 + 11a + ll/3)m2 - 23a/?s] ^ , 



da 



dp 



P ' 
[(a + P)ml - aPs\ml [96(1 - a - /3)2(5 + a + ^) 



a/3 



54(1 - a - /?)(3 + a + /?) 5(1 - a - ^)2(29 + 13a + 13^) 



a 



ap 



90(1 + a + /?) 



[(a + /3)TOc ~ ^1^^ 
7^ 



288(1 - a - ^)a + 60(1 - a - pf + 120a/3 



32(1 - a - Pf[b{l + 2a + 2/3)to2 _ I8a/3s]77i2 
<yGq)f^^-^ ^ "tq(gg^gG'g)(s - lOmg) ^y ,^ Img/,? ^^9scGq) j" 



487r4 



3 X 2^714 



da 



d/3 



(29toc - 9m,)[3(a + P)ml ~ 5a/3s] 6mc(l - a - ^)[3(a + P)ml - 5aPs 



P 

{1 - a - P)[63{a + P)ml ~ 129aPs] 12mqml 



aP 



367r2 



20m; - 2s 



3mqmc{8m'^ — 3s) 
4to? — s 



(A14) 



n 



5 '^'^^^ 967r2 

For the interpolating current Jg^ 

pr*(s) 



2 2(8a2 - a - 20)m2 16(3 - 5a)m 
16)717^ — - 



a(l — a) 



a2(l -a)Af| 



pr^(^) 



{glGG) 

32 X 2i27r6 



da 



/3„ 



d/3 



[(a + P)m\ - aPs]ml 
aP 



48(l-a-/3)2(5 + a + /3) 



a^ 



18(1 - a - /?)(3 + a + /?) (1 - a - /3)2(29 + 13a + 13/3) 



a 



a/3 



+ 18(l + a + /?) 



12[(1 - a - /3)2 + 2a/3][(a + /3)m2 - a/3s]s 16(1 - a - /3)2[5(1 + 2a + 2/3)m2 - 18a/3s]r 



a/3 



mq{qgs<jGq){s - 10m2) 



■\/l — 4m2/s 



{qgsCrGq) f" 



da 



dp 



967r4 ^ 3 X 2^714 

r (13mc + 3TOg)[3(a + P)ml - 5a/3s] (1 - a - /J)(6/3 - 9a)[3(a + P)ml - 5a;3s]m, 



/3 

{qgsaGq){qq ) 
967r2 



a/?^ 

daU8a + ll)A4--2(4"'-^«-^)"^^ 8(3 - Sa^ 



(A15) 



a(l-a) a2(l-a)A/| 

From these results the expressions for the currents Ja^t, J4^, Jy^ and Jg^ can then be obtained conveniently by the 
replacement rric — >■ —ttIc'. 



nic—^ — mc 



> P3(s),P2(s) ^ ^ P4(s),P5(s) 



> P7(s),P6(s) ^ P8(s) 



(A16) 
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4. The spectral densities for the currents with J^*^ — 1+ 
For the interpolating current Ji^: 

,{1 - a - p)[{a + p)ml - aps]^ 



Pi 



{qg^aGq) 



is) 



1 



da / o?/3- 



87r4 



da I dp 



+ 12mcTOg(l - a - /3)[(a + /3)m^ - 3a/3s] 



[(a + /^)m2 - aPs] ( mdl - a- f3)[3{a + P)ml - la^s 



aP 



P 



+ m,[(4 + a + P)ml - 3aPs] 



(GG). X _ {glGG) 

Pi ~ o2 



32 X 2i27r6 



da 



d/3 



[l-a-pfm^ 



96[(a + ;3)m2 - 2a/3s] 5(5 + a + 7/3) [(a + P)ml - aP& 



4(l-a-/3)(5m2 + 2a/3s) 



{qg.aGq) 



5m2(l + a + /3) 



6[(a + /3)m2 — a/Js 

^-^ Vl-4m2/,s-^-^y^^^^^ day^^^^^^ d/3 

77ic(l -a-/^)[27(a + /^)77i2 -49a/?s] Ylm^ml (19771^ - 5mg)[3(a + ^)m2 - 5a/3s] 



/3„ 



/32 



/3 



mc{qq) 



12.2 (4-^ + 2777,-^-^)^1-4^2/., 



(9ff.aGg>{g9>^,,,2^ _ {qg s^jGq) {qq) ,r^^^.j2 , 2(6 - a)r77g ^ 48m 



32x257r2 7o s ' ' ^2^^! 



For the interpolating current J2^: 



(A17) 



pr\^) - ^T^p'n^i pT'i^) = opr(^), pT' (^) = o^r' is) 



pr\s) 



pert I 



ill) I 



Amy 



32 X 2i27r6 



(1-0-^)27772 



48[(a + P)ml - 2a/3s] (5 + a + 7/3)[(a + ^)r772 - aPs 



6[(a + /?)?772 — aPs 



aP 

'mlmq{qgsaGq) 
167r4 



4(1 - a - P){ml - 2aPs) 



a 



P' 

ml{l + a + P) 



,{1 - a - P)[9{a + P)ml - 19aPs] UmqUil {llm.c - mq)[i{a + P)ml - baPs 



P' 



P 



n 



{qgsi7Gq){qq) 



(Af 



{qgs(jGq){qq ) 

32 X 2^712 



dai 3a Mi 



2(6-7a)7772 24777;! 
a{l — a) 



a^Ml 



(A18) 



19 



For the interpolating current Js^: 
1 



(GG)/ N 



3 X 2%6 



da 



d/3- 



(1 - a - /?)[(« + /3)to2 _ Q,/3s]2 J 9(1 + a + /3)[(a + ^)m2 - al3sf 



a 



2«3 



4m2(l -a~ /?)(5 + a + /?)[(q; + /3)ml - 



108mcmq(l -a- P)[{a + P)ml - 3a/3s] ^ , 



da 



dp 



_m r 

87r4 

[{a + ^)m2 - aPs] f 3mc(l - a - /3)[3(a + l3)ml - la/Ss] 
aP 

32 X 2i27r6 



da 



dp 



[{a + ;3)m^ — a/3s]m^ 



5mJ(2 - a - /3)to^ + a/3s] 



a/3 



96(l-a-/3)2(5 + a + /3) 



54(1 -a- /3)(3 + a + /3) 5(1 - a - /3)2(29 + 13a + 13/3) 



+ 90(l + a + /3) 



[{a + P)ml -a/3s]s 



288(1 - a - /3)a + 60(1 - a - pf + 120a/3 



32(1 - a - pf[{h - 8a - 8/3)m2 + 18a/3s]TO^ 



mq{qgaaGq){s + 26m^ 



{qgsCrGq) 



da 



dp 



487r4 " ^ 3 X 277r4 

(1 - a - /3)[87(a + P)ml - 161a/3s]TOc (14 - 71a - 14;3)[3(a + P)ml - 5aPs]mc 



/?2 

9toJ(8 - 3a - 3P)ml + baPs 



aP 



{qq? 
367r2 



o QmamdSmf. — 3s) 

52m2 + 2s ^ " ^ 

4mt — s 



v/1 -4m2/s, 



(A19) 



n 



{qg,,^Gq){qq)fJ. ,J2^ ^ {q9s<yGq} {qq) 
For the interpolating current Jg^ 

(g2GG) 



da<^ (7a + 16)A/^ 



2 2(8a2 - 25a + 4)to2 16(3 -a)r 



a(l — a) 



a2(l-a)Af| 



pr\s) 



32 X 2i27r6 



da 



/3„ 



d/3 



[(a + /3)m2 - aPs]ml 
aP 



48(l-a-/3)2(5 + a + /3) 



+ 18(l + a + /3) 



pf-«^>(s) = 



18(1 - a - /3)(3 + a + P) ^ {1 - a - P)^(29 + 13a + 13/3) 

a a/3 
12[(1 - a - /3)2 + 2a/3][(a + P)ml - a/3s]s 16(1 - a - /3)2[(5 - 8a - 8/3)m2 + 18a/3s]m2 

a/3 a'^ 
^,(,:,..G,Hs + 26mg) ^^^-^_ 1^ £-,^j'-,p 

(1 - a - /3)[3(a + /3)m2 - 13a/3s]mc [(2 + 31a - 2/3)mc + 3am5][3(a + /3)m2 - 5a/3s] 



/32 



n, 



{qgs<yGq)(qq) 



{Ml) 



{qgs(TGq){qq ) 
967r2 



da<^ (8a + ll)A/B- 



a/3 

2 2(4a2 - 9a - 4)to2 8(3 - a)m^ 



a(l — a) 



a2(l -a)M| 



From these results the expressions for the currents Js^, J^fj,, Jjfj, and Js^ can then be obtained conveniently by the 
replacement nic — !• —mc'- 



pi(s) 



> P3(s),P2(s) — ^ p4(s),P5(s) ~ ^ p7(s),P6(s) 



> P8(s) . 



(A21) 



